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PREFACE 


This  Memorandum  presents  an  algorithm  for  choosing  locations  to 
place  air  or  ground  forces  in  order  to  prevent  an  opposing  force  from 
proceeding  through  a  transportation  or  supply  network.  It  is  a  part 
of  continuing  RAND  research  on  the  effectiveness  of  interdiction  and 
is  particularly  applicable  to  infiltration  and  counterinsurgency. 

The  model  has  been  programmed  for  use  on  RAND's  computer,  and 
should  also  be  useful  to  other  organizations,  including  agencies 
interested  in  targeting  strikes  against  lightly  traveled  line-of- 
communications  networks  . 


SUMMARY 


This  Memorandum  presents  an  algorithm  for  determining  where  to 
place  forces  in  order  to  maximize  the  probability  of  preventing  an 
opposing  force  from  proceeding  from  one  particular  node  in  a  network 
to  another. 

The  usual  gaming  assumptions  are  invoked  in  this  model;  namely, 
that  the  strategy  for  placing  forces  is  known  to  the  opponent  and  that 
he  will  choose  a  path  through  the  network  which,  based  on  this  knowledge, 
maximizes  his  probability  of  successful  traverse.  As  given  quantities, 
the  model  requires  a  list  of  the  arcs  and  nodes  of  the  network,  the 
number  of  forces  available  to  stop  the  opposing  force,  and  the  proba¬ 
bilities  for  stopping  the  opposition  at  the  arcs  and  nodes  as  functions 
of  the  number  of  forces  placed  there.  From  this  data,  the  model  cal¬ 
culates  the  probabilities  for  placing  the  force  at  the  arcs  and  nodes 
when  one  force  is  available,  and  the  expected  numbers  ol  forces  to 
place  at  the  arcs  and  nodes  when  multiple  forces  are  available. 

A  computer  program  for  t lie  model  has  been  written  in  Fortran  IV. 
Though  originally  intended  lor  the  1 HM  7044,  it  may  he  adapted  to  lit 
other  computers.  The  program  presently  handles  problems  with  up  to 
300  arcs,  150  nodes,  and  25  foices.  I$y  changing  the  dimension  state¬ 
ments,  these  quantities  can  be  increased  lor  larger  computers,  or  their 
proportions  can  be  modified  for  computers  of  the  same  size. 
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I .  INTRODUCTION 

For  many  years  RAND  has  been  interested  in  problems  of  inter¬ 
diction  and  rnfi It  ration .  Interdiction  has  been  dealt  with  extensively 
from  the  standpoint  of  reducing  the  throughput  capacity  of  a  supply 
network.  The  most  recent  study  owes  most  of  its  theoretical  develop¬ 
ment  to  Wollmer  (6),  operating  models  to  Durbin  (1),  and  operational 
studies  to  J.  U.  Higgins.  This  approach  is  inadequate,  however,  in 
applying  interdiction  to  infiltration  and  counterinsurgency.  The 
reason  is  that  requirements  are  so  small  under  those  circumstances 
that  air  and  ground  bombardment  could  never  hope  to  reduce  a  network's 
throughput  capacity  below  the  level  necessary  to  meet  minimum  require¬ 
ments.  For  these  situations,  interdiction  must  therefore  be  approached 
from  the  standpoint  of  direct  assault. 

In  the  situation  depicted  here,  an  infiltrator  attempts  to  proceed 
from  one  point  to  another  in  a  network.  An  interceptor,  who  may 
possess  one  or  more  forces,  attempts  to  stop  him  by  placing  forces 
along  arcs  or  nodes  that  he  expects  the  infiltrator  to  travel.  His 
problem  is  to  place  his  forces  so  as  to  minimize  the  infiltrator's 
probability  of  successful  traverse.  The  Infiltrator's  problem  is  of 
course  to  select  a  path  that  will  maximize  this  probability.  Both 
problems  can  be  represented  by  a  zero-sum  two-person  game.  As  the 
next  section  will  show,  however,  the  game  matrix  is  very  large  and 
difficult  to  generate.  By  using  an  incremental  approach,  a  solution 
may  be  obtained  much  more  easily  for  the  interceptor. 

This  report  presents  a  model  based  on  sucli  an  approach.  As 
inputs,  it  requires  a  list  of  the  arcs  and  nodes  of  the  network,  the 
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number  of  forces  available  to  stop  the  opposing  force,  and  the  proba¬ 
bilities  for  stopping  the  opposition  at  the  arcs  and  nodes  as  functions 
of  the  number  of  forces  placed  there.  From  this  information,  the  model 
calculates  the  expected  number  of  forces  to  place  at  each  arc  and  each 
node.  For  the  case  where  only  one  intercepting  force  is  available, 
the  expected  values  are  probabilities  of  force  placement;  for  more 
than  one  force,  expectations  are  to  be  interpreted  in  the  most  obvious 
way.  Specifically,  if  the  expected  number  of  forces  to  place  at  a 
particular  location  is  4k,  one  would  place  four  forces  there  always 
and  a  fifth  force  one- fourth  of  the  time. 

The  solutions  obtained  by  this  model  are  always  optimal  in  the 
gaming  sense  (i.e.,  the  infiltrator's  best  chance  of  successful  traverse 
is  as  small  as  possible)  if  the  interceptor  has  one  force  available, 
and  are  optimal  or  nearly  optimal  if  he  has  more  than  one. 
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II.  INTERCEPTOR'S  PROBLEM 

The  network  is  characterized  by  sets  of  elements  called  arcs  and 
nodes.  Nodes  are  points  or  junctions,  and  arcs  are  line  segments 
joining  nodes.  An  arc  that  joins  node  i  to  node  j  is  designated  by 
the  symbol  (i,j).  The  node  from  which  the  infiltrator  starts  is  called 
the  source;  the  one  he  attempts  to  reach  is  called  the  sink. 

It  is  assumed  that  the  infiltrator  and  interceptor  are  each  composed 
of  a  single  force.  Later  on  this  restriction  will  be  relaxed  for  the 
interceptor.  As  is  customary  in  game  theory,  it  is  assumed  that  the 
interceptor's  strategy  will  become  known  to  the  infiltrator.  However, 
this  does  not  necessarily  mean  that  the  infiltrator  will  kno.  Mie 
location  of  the  force.  Specifically,  if  the  interceptor  decides  to 
adopt  a  mixed  strategy  such  as  placing  his  force  at  node  a  with  proba¬ 
bility  ^  and  at  node  b  with  probability  the  infiltrator  will  know 
the  probabilities  but  will  not  know  when  the  force  will  appear  at  a 
and  when  it  will  appear  at  b.  The  infiltrator  will  react  by  traveling 
only  on  source-to- sink  paths  that  maximize  his  probability  of  success¬ 
ful  traverse.  The  interceptor  tries  to  choose  a  location  for  his  force 
that  will  minimize  that  probability.  He  bases  his  choice  on  the 
topology  of  the  network  and  upon  the  vulnerability  of  the  arcs  and 
nodes . 

Note  that  mixed  strategies  are  often  necessary.  To  see  this 
consider  a  network  in  which  the  source  and  sink  are  relatively  invul¬ 
nerable  to  attack  and  any  other  particular  arc  or  node  cjn  be  bypassed. 

In  this  situation,  placing  a  force  at  any  arc  or  node  othei  than  the 
source  or  sink  with  probability  one  would  be  ineffective  since  the 
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infiltratcr  would  then  choose  a  path  that  did  not  Include  the  arc  or 
node  where  the  force  was  placed.  Since  placing  the  force  at  the  source 
or  sink  would  also  be  ineffective,  no  pure  strategy  would  be  desirable. 
To  develop  a  good  mixed  strategy,  it  is  necessary  to  define  quan- 
tities  measuring  the  effectiveness  of  placing  a  force  at  the  various 
arcs  and  nodes,  and  quantities  expressing  the  interceptor's  strategy. 
These  are  as  follows: 

p(i)  =  probability  that  the  infiltrator  will  be  stopped  at 
node  i,  given  that  he  attempts  to  cross  it  and  the 
interceptor  chooses  to  place  his  force  there 

p(i.j)  *  probability  that  the  infiltrator  will  be  stopped  at 
arc  (i,j)  given  that  he  attempts  to  cross  it  and  the 
interceptor  chooses  to  place  his  force  there 

it ( i)  ■  probability  that  the  interceptor  places  his  force  at 
node  i 

Ti ( i , j )  *  probability  that  the  interceptor  places  his  force  at 
arc  (i,j) 

The  quantities  p(i)  and  p(i,j)  may  themselves  include  many 
elements  and  factors.  For  example,  suppose  it  is  determined  that  if 
node  i  is  on  the  infiltrator's  chosen  path  and  the  interceptor  chooses 
to  place  his  force  there,  there  is  a  probability  of  0.8  that  the  inter¬ 
ceptor  will  arrive  there  before  the  infiltrator.  Suppose  further  that 
if  both  were  at  node  i,  the  interceptor's  probability  of  detecting 
the  infiltrator  is  0.9,  and  if  the  two  engaged  in  direct  combat  the 
probability  of  the  interceptor's  winning  is  0.7.  Then  p(i)  would  be 
the  product  of  these  three  factors--0 .8,  0.9,  3nd  0.7 — or  .504. 
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The  probability  that  the  infiltrator  can  successfully  cross  a 
particular  node,  a^,  is  1  -  tT(a^)p(a^),  and  that  of  crossing  a  partic¬ 
ular  arc,  (a^a^)  is  1  -  n(a^,a  )p(a^,aj)  .  Assuming  these  probabilities 
★ 

are  independent  and  the  infiltrator  attempts  to  reach  the  sink  by  the 
path  a^,  a^,  his  probability  of  successful  traverse  is 

n  n 

(1)  K  *  n  [  1  -  rr(a^)  x  pCa^]  x  H  [l  -  "(at»ai+1)  x  P(ai»ai+1)l* 
i«i  i=l 


This  quantity,  K,  will  be  referred  to  as  the  value  of  the  path  a^  , 

a  .  Given  the  p(i),  p(i,j),  rr(i),  and  rr(i,j),  the  infiltrator's 
n 

problem  is  to  find  a  source-sink  path  of  maximum  value.  The  interceptor's 
problem,  given  the  p(i)  and  p(i,j),  is  to  choose  tt ( i)  *•  0  and  n(i,j)  £  0 
such  that  the  maximum  value  of  all  source-sink  paths  is  minimized, 
subject  to  the  constraint  I  n(i)  +  I  tt ( i ,  j )  s  1. 

The  problem  of  finding  the  rr(i)  and  n(i,j)  is  really  one  in  game 
theory.  It  can  be  represented  by  a  game  matrix  and  solved  by  linear 
programming.  The  linear  program  solution  would  include  an  optimal 
strategy  for  the  infiltrator  as  well  as  the  interceptor.  The  value  of 
the  game  would  be  equal  to  both  the  infiltrator's  maximum  guaranteed 
probability  of  successful  traverse,  and  one  minus  the  interceptor's 
guaranteed  probability  of  stopping  the  infiltrator.  However,  the  game 
matrix  would  require  a  column  for  every  possible  source-sink  path  and 
a  row  for  each  of  the  Interceptor's  pure  strategies.  (The  number  of 
these  pure  strategies  becomes  exceedingly  large  when  the  interceptor 
is  allowed  multiple  forces.)  Hence,  this  formulation  is  impractical. 


While  this  assumption  is  not  strictly  true,  it  will  be  shown 
later  that  this  leads  to  no  inaccuracies  for  the  one-intercept i ng- force 
case . 
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The  incremental  approach  to  be  presented  avoids  these  problems 
of  enumeration  and  the  resulting  large  number  of  variables  for  the 
linear  program.  Furthermore,  it  yields  an  optimal  strategy  when  the 
interceptor  has  but  one  force,  and  an  approximately  optimal  strategy 
when  he  lias  more  than  one  force. 
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III,  AN  INCREMENTAL  APPROACH  FOR  ONE  INTERCEPTING  FORCE 

This  section  solves  the  interceptor's  problem  by  a  marginal 
analysis  or  st eepest-ascent  approach  when  the  interceptor  has  hut  one 
force . 

Specifically,  the  approach  is  as  follows.  Initially,  all  n(i)  and 
rt ( i ,  j )  are  assigned  values  of  zero.  Then  they  are  increased  by  small 
amounts  An(i)  ^  0  and  An(i,j)  -  0.  The  proportions  for  the  An(i)  and 
Arr(i,j)  are  such  that  the  additional  force  allocation,  E  An(i)  +  E  An(i,j 
divided  by  the  decrease  in  maximum  path  value,  is  minimized  as  these 
two  quantities  tend  toward  zero.  This  is  equivalent  to  maximizing  the 
decrease  in  maximum  path  value  per  unit  of  additional  force  allocation. 
The  rt(i)  and  tt ( i ,  j )  are  then  increased  again  in  the  same  manner  until 
they  sum  up  to  one. 

For  a  specific  path,  a^,  ...,  a^,  increasing  ufa^)  by  unfa^)  replace 
the  factor  [l  -  n(a^)p(a^)j  in  Fq  .  (1)  by  [l  -  (n(a^)  +An(a  ))p(a  )]  and 
reduces  the  value  of  the  path  from  K  to  K  -  uK,  where 


K 


AK 


K 


Ln(ai) 


+  An  ( 

n(ai  ) 


a . 
t 


x 


)J  X  l p(a^) J  j 
Top  )• 


Solving  for  the  quotient  of  the  additional  iorce  allocation  and  the 
decrease  in  path  value,  the  following  expression  is  obtained: 


If  force  allocation  were  increased  at  an  arc  instead  of  a  node, 


^nUd  .  l  r  i 

AK  K  |^p(a.) 


the  same  argument  would  yield  an  expression  lor  the  additional  lorce 
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allocation  divided  by  decrease  in  path  value  identical  to  the  preceding 
one,  with  arc  probabilities  substituted  for  node  probabilities.  Speci¬ 
fically,  for  arc  this  wou^  be 


An(a. 


,ai+l) 


AK 


*[' 


P(ai»ai+1) 


-  n(a. 


'W 


Note  that  in  order  to  obtain  any  decrease  in  maximum  path  value, 
the  An(t)  and  An(i,j)  must  be  strictly  greater  than  zero  along  a  subset 
C,  of  the  arcs  and  nodes  that  intersects  all  paths  of  maximum  value. 

An  efficient  allocation  would  require  that  all  members  of  C  remain  on 
maximum  value  paths,  for  otherwise,  lessening  the  force  increase  at 
an  arc  or  node  not  on  a  maximum  value  path  by  a  small  amount  would 
decrease  XAn(i)  +  Z&rr(i,j),  while  not  affecting  the  maximum  path  value. 
Then,  of  course,  this  small  amount  of  force  could  be  redistributed 
among  arcs  and  nodes  on  maximum  value  paths  to  get  a  strict  improvement 
Thus,  the  An(i)  and  Att(1,J)  assigned  to  the  chosen  C  must  be  such  that 
the  AK's  are  equal  for  each  of  its  members.  In  other  words,  they  must 
be  proportional  to 


(2) 


a) 


pU) 


-  T7  (  i) 


b) 


1 


p(i. j) 


rr  •  TT(i*J) 


The  constant  — ,  where  K  is  the  maximum  path  value,  was  dropped  in 
Expression  (2),  of  course.  The  total  increase  in  force  allocation 
per  unit  decrease  in  the  value  of  a  maximum  probability  path  would  be: 
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(3>  kIIju) -"<i>]  +  k  I  bubr •  n(t-J)] • 

t«C  (i,j)eC 

Thus,  the  problem  is  reduced  to  finding  a  subset  of  arcs  and  nodes 

intersecting  all  maximum  probability  paths  that  minimizes  Expression  (2). 

Initially,  all  paths  have  value  one  since  all  rr ( i )  and  tt ( i ,  j )  are 

zero.  Thus  C  is  required  to  be  a  set  blocking  all  source-sink  paths. 

If  the  nodes  and  arcs  are  assigned  capacities  equal  to  the  quantities 

of  Expression  (2a)  and  (2b),  the  problem  of  finding  C  is  reduced  to 

★ 

one  of  finding  a  minimum  cut.  However,  this  is  equal  to  the  value  of 
the  maximum  source-sink  flow  and  can  be  found  by  the  maximum  flow 
algorithm  of  Ref.  3. 

Note  that  once  additional  force  is  allocated  to  the  minimum  cut, 

C,  the  capacities  of  all  arcs  and  nodes  in  C  decrease  by  the  force 
Increase  while  the  capacities  of  all  other  arcs  and  nodes  remain  the 
same.  Thus,  the  value  of  C  will  be  reduced  by  L2m(i)  +  Z&Tr(i,j)  while 
all  other  cuts  will  decrease  by  amounts  that  do  not  exceed  this.  Hence, 

C  remains  minimum  and  the  entire  unit  of  force  may  be  allocated  to  C. 

The  tt ( i)  and  n(i,j)  may  therefore  be  solved  for  as  follows: 

1.  Assign  all  nodes  capacities  of  l/p(i)  and  all  arcs  capacities 
of  l/p(i,j)  (since  all  tt ( i )  and  rr ( i , j )  are  zero). 

2.  Maximize  flow  from  source  to  sink. 

3.  Let  C  be  the  minimum  cut  set  and  V  its  value.  Set  n(i)  = 
p(i)/V  for  i  €  C  and  r(i,j)  =  p(i,j)/V  for  (i,j)  e  C.  Set 

★ 

Essentially,  a  cut  set  is  a  set  of  arcs  and  nodes  blocking  all 
source-sink  paths.  Its  value  is  the  sum  of  the  capacities  ol  its  arcs 
and  nodes,  and  the  minimum  value  of  all  cuts  is  equal  to  the  maximum 
value  of  the  source-sink  flow. 
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all  other  n(i)  and  tt  (  i ,  j )  equal  to  zero. 

Note  that  since  only  one  cut  is  obtained,  each  path  of  maximum 
value  has  force  allocated  to  only  one  of  its  arcs  or  nodes,  and  both 
its  value  and  traverse  probability  are  equal  to  the  probability 
associated  with  that  particular  arc  or  node.  Hence,  the  independence 
assumption  does  not  lead  to  inaccuracies  for  the  one-intercepting- 
force  case.  Finally,  the  maximum  path  value  is  1  -  ~,  where  V  is  the 
value  of  C  found  in  step  3. 
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IV.  MULTIPLE  INTERCEPTING  FORCES 

The  incremental  procedure  of  the  last  section  may  be  modified  to 
handle  the  more  general  case  in  which  there  is  more  than  one  intercepting 
force.  The  quantities  measuring  tiie  effectiveness  of  placing  forces 
at  the  various  arcs  and  nodes  and  those  defining  the  interceptor's 
strategy  must  be  redefined  as  follows: 

p(i)^  =  probability  that  the  infiltrator  will  he  stopped  at 

node  i,  given  that  he  attempts  io  cross  it  and  the 
interceptor  lias  placed  k  of  his  lorces  there; 

p(i,j)  =  probability  that  the  infiltrator  will  be  stopped  at 
arc  (i,j),  given  that  he  attempts  to  cross  it  and 
the  interceptor  places  k  of  his  forces  there; 

rt(i)  =  expected  number  of  forces  the  interceptor  places  at 
node  i ; 

n(i,j)  =  expected  number  of  lorces  the  interceptor  places  at 
arc  (i ,  j)  , 

where  P(Oq  and  p(i,j)g  are  both  identically  zero.  Also  let 

n(i)  =  tt ( i )  -  [n(t)  ] 
n(i,j)  =  n(i , j )  -  [n(i,j)] 

Note  that  n(i)  and  fr(i,j)  are  merely  the  fractional  parts  ol  n(i)  ami 
rr(i,j).  It  will  be  assumed  that 

p<1)k+i  ’  p(1)k  p(i)k  • 

P(i.j)k+1  "  P(i.j)k  '  P(*.j)k  -  P(i.J)k_r  k  '  1- 


Mathematically,  this  is  essentially  a  convexity  assumption  on  the  value 
of  forces  at  the  arcs  and  nodes.  Physically,  it  expresses  the  lact 
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that  multiple  forces  may  result  in  overkill.  This  assumption  assures 

that  for  a  given  value  of  rr(i),  the  probability  associated  with  node  1 

Is  minimized  by  allocating  [rr(i)]  +  1  forces  at  i  with  probability 

n(i),  and  [n(i)]  forces  with  probability  1  -  n(i) .  In  other  words,  if 

n(i)  =  3  2 ,  n(i)  would  be  realized  by  assigning  three  forces  at  i  half 

the  time  and  four  half  the  time,  as  opposed  to  such  a  policy  as  allocating 

two  forces  half  the  time  and  five  half  the  time.  A  similar  result  holds 
★ 

for  the  arcs.  Thus,  only  policies  of  this  type  need  be  considered 
and  the  n(i)  and  n(i,j)  completely  define  the  interceptor's  strategy. 

As  in  Eq .  (1)  of  Sec.  II,  the  value  of  a  path  is  still  the 
product  of  the  probabilities  associated  with  its  nodes  and  arcsj 
however,  the  probability  associated  with  node  i  is  now 

S(U(i  *  *  (1  '  S<t))(1  -  P<1>ln<l>]> 

or 

1  '  P<UU<L>]  '  S<1)(P(1)[„(I)]+1  -  P(1>ln(l>]) 
and  that  with  arc  (i,j)  is 

If  a^,  ...,  a^  is  a  path  of  value  K,  and  n(a^)  is  increased  by 

An(a^),  reducing  the  path  value  to  K  -  AK,  the  expressions  for  K  and 
K  -  AK,  as  before,  differ  only  in  the  factor  for  node  a^.  Thus,  for 
sufficiently  small  An(a^),  the  expression  for  the  quotient  ol  the 
force  increase  and  path  decrease  is 

★ 

See  Appendix  C . 
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ATr(at) 

~AK 


J, 

K 


1  '  p(V[tt(V3 _ 

P(*i,[n(a1)]+l  ‘  p(at)[n(a1)]  "<*i) 


and  decreasing  path  value  by  increasing  force  at  an  arc  instead  of  a 
node  yields  the  analogous  expression, 


An(at>a  )  l 
AK  K 


1  ~  P(afaj)[n(ai,a.)] 


P^al’aj^[n(ai,a^)]+1  "  p^aL ’ 3 j J [n(ai ,a. ) j 


-  "(a^a^) 


As  in  Sec.  Ill,  force  must  be  allocated  along  all  arcs  of  some 
cut  set,  and  the  problem  is  reduced  to  finding  a  minimum  cut  set,  which 
in  turn  reduces  to  one  of  finding  a  maximum  source-sink  flow.  The 
node  and  arc  capacities  are  now 


(A) 


a)  c(i) 


LlflVlnUJl-  ; _ 

P<l)[nU)]+l  '  p(l)ltT(i)] 


n(i) 


b)  c(i.j) 


1  •  '<1'J)1.(U)1 
p<i,))[n(t,j)]+l  '  p(l,j)ln(i,j)] 


n(i,  j)  . 


Note  that  the  c(i)  and  c ( i , j )  are  no  longer  decreasing  functions 
of  tt ( i )  and  rr(i,j)  except  over  intervals  whose  endpoints  possess  the 
same  integer  part.  Specifically,  c(i)  is  a  decreasing  function  of 
tt(1)  for  0  s  rr(i)  <  1 ,  1  s  n(i)  *-  2,  etc.  However,  c(i)  increases  at 
the  points  n(i)  **  1,  2,  etc.  Thus,  when  a  minimum  cut  set  is  found  ami 
additional  force  is  allocated  along  its  arcs  and  nodes,  that  cut  remains 
minimum  provided  the  forces  on  these  arts  and  nodes  do  not  leach  or 
exceed  their  next  highest  integer  values.  This  requires  that  t  lie 
Incremental  approach  be  modified  as  lollows. 
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Starting  with  all  n(i)  and  rr(i,j)  equal  to  zero,  capacities  are 
assigned  to  the  arcs  and  nodes  o£  the  network  which  are  equal  to  the 
quantities  defined  by  Expression  (3).  Flow  is  maximized  from  source  to 
sink  and  the  n(i)  and  n(i,j)  corresponding  to  the  arcs  and  nodes  of  the 
minimum  cut  set  are  increased  by  amounts  An(i)  and  An(i,j)  which  are 
proportional  to  their  capacities.  The  constant  of  proportionality,  M, 
is  the  smallest  possible  such  constant  that  will  either  increase  some 
rr(i)  or  rr  ( i ,  j )  to  the  next  highest  integer  or  will  increase  Lrt(i)  + 
£n(i,j)  to  n.  If  the  latter  does  noL  happen,  the  new  values  of  n(i)  and 
r(i,j)  are  used  to  calculate  new  capacities  and  flow  is  maximized  again 
to  obtain  a  new  cut.  The  process  is  repeated  until  £n(i)  +  £rr(i,j)  *  n. 

Specifically,  the  algorithm  for  the  general  case  of  n  intercepting 
forces  is  as  follows: 

1.  Set  all  n(i)  *  0  and  ail  rt(i,j)  =  0. 

2.  Assign  the  arcs  and  nodes  of  the  network  capacities  of  c(i) 
and  c(i,j)  respectively,  where  c(i)  and  c(i,j)  are  as  defined 
by  Eq .  (4) . 

3.  Maximize  flow  from  source  to  sink  and  let  C  be  the  arcs  and 


nodes  of  the  minimum  cut  set  found. 

4.  Compute  M.  =  min  ^ — "  7! ^ 

1  i€  C  t(L) 


M,  -  min  U  ~  "U,. 
2  (l.j)ec 


M.  =  n  -  i-n(i)  -  In(i,  j) _ 

L  c(i)  +  L  c  ( i  ,  j ) 
ieC  (i,j)eC 


Let  M  = 


min  (M1,M2,M  ).  Set  An(i)  =  Mc(i)  if  ieC  and 


Att(1,  j)  =  Mc(i,  j)  if  (i,j)eC.  Set  all  other  An(i)  and 
( i» j )  equal  to  zero. 
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6.  Set  n(i)  *  n(l)  +  An(i)  and  rr(l,j)  =  rr(i,j)  +  Art(i,j)  for  all 
i  and  (t,j).  If  £  rr(i)  +  £  tt ( l , j )  <  n,  go  back  to  step  2. 
Otherwise  terminate. 

In  step  4,  is  that  constant  of  proportionality  which  would 
increase  some  rr(i)  to  an  integer,  the  one  which  would  increase  some 
tt ( i ,  j )  to  an  integer,  and  the  one  which  would  increase  £  n(i)  +  £  -n ( i ,  j ) 
to  n . 

At  termination  the  rr ( i )  and  rr(i,j)  represent  the  expected  number 
of  forces  to  place  at  node  i  and  arc  (i,j),  respectively.  The  policy 
will  be  to  place  [n(i)]  +  1  forces  at  node  i  witli  probability  n(i)  and 
[n(i)]  forces  with  probability  1  -  n(i).  The  value  of  the  resulting 
maximum  probability  path  can  be  calculated  by  the  algorithm  of  Appendix  D. 

If  the  iterations  yield  a  sequence  of  but  one  minimum  cut  set 

(i.e.,  if  steps  2  through  6  are  performed  only  once),  the  solutions 

found  are  optimal,  as  in  the  case  of  one  intercepting  force,  and  also 

the  independence  assumption  on  the  arc  and  node  probabilities  leads 

to  no  inaccuracies  in  maximum  path  value.  However,  if  the  procedure 

yields  a  sequence  of  several  cuts,  neither  optimality  nor  independence 

can  be  guaranteed.  Nevertheless,  the  nature  ol  t he  strategy  indicates 

that  independence  can  be  violated  only  by  the  fractional  parts  of  the 

n(i)  and  tt ( i , j )  .  The  maximum  path  value  can  be  kept  track  of  during 

the  course  of  the  algorithm.  Specifically,  if  K  is  the  maximum  path 

value  at  the  beginning  of  an  iteration,  then  the  maximum  path  value 

at  the  end  of  the  iteration  is  K ( 1  -  M) ,  where  M  is  that  found  in  Step 
★ 

5.  Of  course,  K  =  1  at  the  beginning  of  the  first  iteration. 

* 

See  Appendix  D. 
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Appendlx  A 

COMPUTER  PROGRAM  (FLOG) * 

A  computer  program  for  the  model  was  written  in  Fortran  IV  for 
use  on  the  IBM  7044.  It  can  easily  be  adapted  for  use  on  other  high 
speed  computers.  The  main  inputs  are  the  total  number  of  forces 
available  to  the  interceptor,  and,  for  each  arc  and  node,  his  probability 
of  stopping  an  infiltrator  attempting  to  cross  it  as  a  function  of  the 
number  of  forces  placed  there.  The  outputs  are  the  expected  number  of 
forces  to  place  at  the  arcs  and  nodes  and  the  infiltrator's  maximum 
probability  of  successful  traverse. 

The  program  presently  handles  problems  with  up  to  300  arcs,  130 
nodes,  and  25  forces.  These  capabilities  can  be  modified  to  meet 
individual  needs  by  changing  the  dimension  statements. 

In  changing  the  dimension  statements,  note  that  the  subscripts 
of  the  Nl,  N2 ,  X,  CA,  PA,  and  DPA  arrays,  ami  the  first  subscript  of  the 
U  array,  are  all  equal  to  the  maximum  number  of  arcs.  The  subscripts 
of  the  NL1,  NL2 ,  Y,  CN ,  PN,  and  DPN  arrays,  and  the  first  subscript 
of  the  V  array,  are  all  equal  to  the  maximum  number  of  nodes.  The 
maximum  number  of  forces  is  equal  to  the  second  subscripts  of  the  U 
and  V  arrays. 

★ 

The  program  of  this  section  was  written  by  Steve  Giaseman  and 
Richard  Clasen. 
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FLOG 

DATA  SUBMISSION  INSTRUCTIONS 


Card  1: 

Col. 

Data 

1-3 

Total  number  of  nodes. 

4-6 

Total  number  of  arcs. 

7-9 

Number  of  source  node. 

10-12 

Number  of  sink  node. 

13-15 

Total  number  of  forces 

involved  . 

16 

Output  flag.  0  *  print 
of  problem. 

output  only  at 

1  =  print  output  after  each  iteration. 

17  Input  flag.  0  *  no  data  cases  following, 

1  =  case  following  present  case. 

II.  Group  1: 

Punch  the  number  of  the  first  node  of  each  arc  as  follows: 

a .  Start  in  col.  1 . 

b.  Three  columns  per  entry. 

c.  Ten  entries  per  card. 

III.  Group  2: 

Punch  the  numbe r  of  the  last  node  of  each  arc  in  the  same 
format  as  group  1  . 

IV.  Group  3: 

Punch  t he  arc  input  probabilities  as  follows: 

a .  Start  in  col  .  1 . 

b.  Six  columns  per  entry. 

c.  Ten  entries  per  card. 

As  a  guide,  consider  the  following:  with  eight  forces  and  ten  arcs, 
group  3  would  consist  of  eight  cards,  each  with  ten  entries.  The  first 
card  would  contain  arc  1  with  1  through  8  forces,  and  arc  2  with  1  and 


2  forces,  etc. 
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V.  Group  4: 

Punch  the  node  input  probabilities  in  t he  same  format  as  group  3 
VI .  Notes 

1.  All  card  entries  are  right- justi fied  . 

2.  Data  deck  must  adhere  to  group  order  in  above  instructions. 

Each  node  must  be  identified  by  a  number  not  to  exceed  three 
digits.  Output  may  be  requested  at  the  end  of  each  iteration  or  maxi¬ 
mum-flow  problem,  or  may  be  requested  only  at  the  conclusion  of  the 
entire  problem.  Multiple  sets  of  data  may  be  run.  If  this  is  done, 
a  0  must  be  placed  in  column  17  of  the  first  card  of  the  last  case, 
and  a  1  in  column  17  of  the  first  card  of  all  other  cases. 

The  output  format  and  the  program  follow. 

OUTPUT  FORMAT 

SOURCE  NET  FLOW  SINK  MAX  PATH  VALUE 


XXX  XX.XXXXX  XXX  x.xxxxx 


NODE 

FLOW 

CUT  SET 

FORCES 

XXX 

X .xxxxx 

X 

x.xxxxx 

ARC 

FLOW 

CUT  SET 

FORCES 

XXX 

X .xxxxx 

X 

X .xxxxx 

A  one  in  the  cut  set  column  indicates  that  that  particular  node 
or  arc  is  in  the  minimum  cut  set  and  hence  will  have  its  number  ol 
forces  increased.  A  zero  indicates  it  is  not  in  the  cut  set.  The 


other  headings  are  self-explanatory. 


1 
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F  IBFTC  FLllf. 

C 

t  TO  EXAMINE  bi.M.RAL  A  TTACKFR  -  F  VADFR  NETWORK. 

C  INPUTS  *  A  T I  ACKER  •  S  AKC  AND  NCDL  FORCE  PLACEMENT  PROB. 

C 

DIMENSION  MII3CM,  N2I30U),  01300,25),  VN50,?5),  XI300). 
lYIISOI,  CM3L0I,  NL1U50),  NL2U50),  0RJI1),  NCSU50),  AtSISO.SO), 
2PNI150),  CNI  150  »  PAI300),  DPNIL50),  0PAI300) 

INTEGER  ACS 
C 

C  INITIALIZATION. 

C 

IH  00  IN  I  =  l, SC 

Y  (  I  »  =  0. 

Nil  (  I)  =  0 
NL  2 (  I  )  =  0 

NCS I  I )  =  (; 

opn (  n  =  o. 
no  ?c  j  -  use 

20  ACS (  I,JI  *  0 
IN  CONTINUE 

DO  21  I  =  I,  ICC 
PAIN  =  0. 

PN (IT  =  0. 

CAIN  =  0. 

21  CNI  I)  =  U. 

no  22  I  =  I,  ISO 
XI N  =  C. 

LPAI N  *  0. 

22  CONTINUE 
PV  -  1. 

C 

C  READ  INPUT  PAPAPF  TEWS. 

C 

READ  lOOu.  M,N»  I  SRC  » I SNK, IF  » IFLAGt  JFLAG 
READ  2000,  INIIN,  I  *  1  ,NI 
READ  2000,  IN2I  N,  I  =  1,NI 
READ  30C0,  I  I U I  I , J  > , J  =  1  ,  I F I  ,  I  =  1  , N » 

READ  3000,  l(VII,J),J  =  l  ,  IF )  ,  I  =  1,P) 

F  =  NOATI  IF  » 

C 

CP*  NCMfUR  OF  NODES. 

C  N  NUMBER  OF  ARCS. 

C  I  SRC  =  SOURCf  NOOE. 

C  ISNK  =  SINK  NODE. 

C  IF  =  number  uf  forces  involved. 

C  IFLAG  =  PRINT  FLAG.  0  =  END  PRINT  ONLY,  I  =  RUNNING  ACCOUNT. 

C  JFLAG  =  INPUT  FLAG.  0  *  1  DATA  SET,  l  *  MULTIPLE  DATA  SETS. 

C  Ml  N  =  ID  nr  FIRST  NOOE  of  ARC  I. 

C  N2(  N  =  ID  OF  LAST  NODE  OF  ARC  I. 

f.  U I  1  »  J  )  -  PROB  Afl  .  of  ARCS 

(.  V(l,J)=  PKObAB.  Of  NODES 

C 

C  C FAN Of  ARC  AND  NODE  PROP.  INTO  CAPACITIES. 

1. 

no  i  j  *  i,n 

I  LAI.II  =  1  ./U(  J,  N 
n  ?  j  =  i ,  p 

?  C.MJT  =  1  .  /V I  J  ,  N 


L 


J 
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C  SEND  COLUMN  HI  ALGORITHM. 

C 

100  CO  101  I  =  l.M 

101  NCS!  I)  *  0 

DO  10?  I  =  l.N 
Nil  =  N  1  (  I  ) 

N2I  =  N?lll 

102  ACSIN1  I  t iN 2 1  I  =  0 

CALL  FLOWMX! v.N,  I  SRC  .  I  SNK.Nl  ,N2  •  C  A  ,  X  #CN  ♦  Y  •  NL  l  ♦  NL2  .  (IF*  J  ) 
C 

C  UPON  RETURN,  FIND  ARCS  AND  NODE S  IN  CUT  SET. 

C 

00  *  I  -  l»M 

3  IFINLK  I).Eu.P.AN0.NL2(  I ) • NC .0 )  NCS(I)  =  1 
CO  4  I  =  l.N 

Nil  =  Nil  I) 

N2I  =  N2(  I  1 

IF(NL?(N1!  ).FC.C.AN0.NLl(N2l  l.NE.O)  ACS(NII,N2I)  =  l 

4  I  F  (  Nl.  2  (  N2  I  1  .  CL'  •  0  .  AND  .NL  l(Nll).Nf-.O)  ACS ( N  l  I  ,N2 I )  =  l 
C 

C  CALCULATE  MOST  LIMITING  ARCS  AND  NODES  IN  CUT  SET. 

C 

SM2  =  99'D. 

CO  ■)  I  =  l.N 
Nil  =  Nil  I) 

N 2  I  =  M?(  I  ) 

IFIAC  SI  NlUN2lUNF.il  GO  TO  b 
C 

AM  2  =  I  1  FLOAT!  IF  I  XI  PAI  I  U.  00002)  1  ♦  1.  )  -  P  A I  I)  )  /C  A  I  I  ) 

IF! AM2.GE.SM2)  GO  TO  b 

SM2  =  AM? 

5  CONTINUE 
C 

C  SM2  IS  NOW  Mf'ST  LIMITING  ARC  IN  CUT  SFT. 

C 

SMI  =  'TSNN . 

CO  6  I  =  l.M 

IFINCSI  D.NC.  1)  GO  TO  6 

AM  1  =  I (FLOAT!  IF  IX! PNl  1 )♦. 0  0002 ) ) ♦ I . ) -PN! I )  ) /CN I  I ) 

IF! AMI. GE. SMI)  GO  TO  6 
SMI  *  AMI 
fc  CONTINUE 
C 

C  SMI  is  NOW  or  ST  LIMITING  NODF  IN  CUT  SFT. 

C 

SUMNP  =  0. 

SUMNC  =  U. 

CO  7  I  =  1,M 

SUMNP  =  SUMNP  ♦  PN(  I  ) 

IF  (  NC  S  (  I  UNI.  I  )  GO  TO  7 
SUMNC  =  SUMNC  ♦  CN(I) 

7  CONTINUE 
SUMAP  =  0. 

SUMAC  =  0. 

DO  8  I  =  l.N 

SUMAP  =  SUMAP  ♦  PA! I  I 

Nil  *  Nl!  I  ) 

N2I  =  N ?!  I  ) 

IF(ACS(NlUN?l).Nfc.l)  GO  TO  8 
SUMAC  -  SUMAC  ♦  CAII) 


-22- 

ft  CONI  IHl>L 
C 

SMI  =  (  J  -  SGMNP  -  SOWAi*)  /(SGMNC  ♦  SOM  AO  ) 

c 

C  l  INI'  M'TSI  LI  III  NO  M  OP  ALL  M*S. 

r 

ALIf  :  AMINlf  SH  •SM2»SMtt 

0 

C  CALCULUI  MAXIMUM  PATH  VALUF. 

(. 

PV  =  PV  ♦  (  l.  -  ALIM) 

t. 

C  INCRFMLNI  i.MA-PI  ARRAYS. 

C 

CO  9  I  -  If' 

I F ( NO  S l  I  )  .Nl  .1)  GO  TO  9 
PNI  l  )  *  Mill  ♦  ICN<  I  )  *  ALIM 
9  C(  NT  !NU( 

co  ic  i  =.  l ,  i 

MI  -  M(  !  ) 

Ml  -  N2I  I) 

IF  (  ACM  NI  1  ,  i.M  )  .NT  .  1  )  GO  TO  10 
PAI  I  )  r  (Mil  I  ♦  I C  A {  I)  ♦  ALIM) 

IC  COM  I Nli! 

C. 

C  ItSI  FIX  ION!  .  IF  NO  f  f  COMPUTE  KF  N  CAPACITIFS.  PRINT  STAII. 

L 

SPNI  =  0. 

no  <>1,1  i  =  i  •  m 

20  1  SPNI  =  hPNI  ♦  PNI  I  ) 

SPA  I  =  0. 

no  2C2  I  =  I  *  N 

20?  SPA  I  =  SPA  |  ♦  PAII  ) 

TFST  =  ,PN|  ♦  SP  A I  ♦  .CCD? 

IF!  TFST.GF.F  .OR. IFLAG.EG.l )  GL’  TC  13 

?oc  no  ii  i  =  i  f  M 

IFINCSI  I  )  .  N*  .  1  )  GO  TO  i  1 
J  =  1 1  I  X ( PNI  I  )  ♦  .0002  ) 

IF (J. I  0.0)  GO  TO  100 

C 

C 

CNII)  =  III. -VII.  J  )  )  /  (V(i»JM)-VI!tJ)))  -  PNt  I  I  ♦FLOAT  (J) 
GC  TO  1  I 

ICO  CM  I  )  =  CNII)  -  PNI  I  I 

11  CONTINOF 

CC  1?  I  =  I  » N 

Nil  =  N  I  (  I  ) 

N  2  I  *  N  2 1  I  ) 

IF  I  AC  SIN)  I.  *|2  I  I.NT.l)  GO  TO  1? 

J  =  iriXIPAl  (I  ♦  . CC02 ) 

IFIJ.M..U)  Gi.  TO  <t  C  0 

( 

CAM)  -  III. -Oil,  J  )  )  /  (  UI  I  #)♦  1  )  -U  I  I  »  J)  )  I  -  PAI  I)  ♦FLOAT  I  J  ) 

GO  T (  12 

ACO  CAM)  =  (.  A I  I  »  -  PAII) 

(. 

I  2  CON  T  I  Mil 

r.u  ii  i oo 


o  r>  o 
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C 

13  w«  1 i r  ( f»,4ccu ) 

w*R  I T  r  (6t«roi )  iw.,  nnj,  isnk,  pv 
WRI  rHr.,4C0.M 
CO  14  I  =  l.M 

14  WKlTt  (6.4CC4)  I,  Y(l),  M  S  (  1)  ,  PM!) 

KR  I  TG  !  b,4  00M 

L  Cl  IS  I  =  |,M 

Is  1 1  -  Mil) 

N 2  1  =  J2I  I  ) 

IF(  X(  l>  .1  T.  :!.)  r.i’  TO  lb 

WHITE  (6,401.  1)  Ml,  N2I,  X(I),  A  l  S  (  N!  I  ,  N  2 !  )  ,  PA  111 

15  ClJNMM.l 
GO  O'  1/ 

lb  XI  =  XII) 

A  VX  =  ALMXI) 

WRlTHb,4CC3)  N2I,  Ml.  AVX,  AG  SI  Ml.',  I),  P  A  (  I  ) 

GO  11)  l» 

17  1M  ICSI.ll. I  )  GO  III  /{IV. 

c 

c  rcm,  get  on-  kagoinf. 
c 

Milt  Ib.SITO) 

IH  JFl  M,.u.  1  )  GO  to  11! 

CALI  CXI! 

FORMA  IS. 

1000  CORMAM  41  1,  I  I, 21  1  1 
2000  FORMAT!  10 1  3) 

3000  FORMAT!  lOKb.M 

4CCC  FORMAT!  lH|,r.  IX.bHSnURCl  ,3X,HHXjF  I  F  I  PW  ,  3  X  ,  4HS  1  NK  ,  24X  ,  l  r»MM  A  <  .  t»A  I M  V 

l ALUF  1 

4001  FORM  A  I  (  IMG,  »4X ,  I  3 ,  X  ,  F  l  2  .  >  .  J  X  ,  1 1 ,  25  X  ,  F  1  2  .  5  /  /  ) 

4002  FORMATIIM  ,42X,  4  HMJO  F  .  IPX, 4  HF  LOW,  I0X,  7HCUT  STT,  7  X  ,  oMH  Ml.  I  5  /  t  I 

4003  FORMATIIM  ,4iX,  I  3  ,  I M  »  ,  13  «  3  X  ,  F  1  2  .  5  ,  1  1  X  ,  I  2  ,  »X  ,  F  1  ?  .  5  ) 

4004  FORMATIIM  ,42X,  13  .  5  X  ,  F  1 2  .  5  ,  1 1  X  ,  1  2  ,  3  X  ,  F  1  2  .  M 

4005  FORMAT!  IhO, //  44X,  30  ARC  ,  l  OX  ,  4MF  l  Uh ,  t  OX  ,  7HC  J  T  SET,  7  X  ,  b*  t  :  *!•  i  S  /  /  I 

5C00  FORMA!  IIHC,  b/t.  140FND  OF  PHOHLE y  I 

END 

I  IHFT(.  FLONMX 

SlIP  ROUT  INI  FI  l,W*X|  NOOE  S,  ARCS,  SRC  ,  SNk  ,  I  ,  J  ,H  I  ,  F  LOW,  NC «  NC  ,  NA  ,  N!i  ,  t'F  I  I 
INIFGFR  NOOf  S.ARGS. SNK, SRC, I  I  2000 1 , J I  2000 ) ,H| 12000) , U OK  I 200U  1 
INTEGER  NCI  I  ICO)  ,  NF  I  ICCOl  ,NA|  IOUOI  ,NH|  I'lOO)  .OBJ 
C 

RCHGIKI  =  FL  MIM/FI 
ACHG!  x.  1  =  l  l  •AMI'llt  X,F?) 
tl  =  2.**1G 
E 2  =  2. **24 
EQUIVALENCE  II  ,P) 

DO  1  Jh  I,  ARCS 
L  =  HU  JJ  ) 

1  H  1 1 J  I  )  =  ACHGtP) 

CO  2  11=  I, NODES 
l  =  NC !  II) 

2  NC!  I  I  )  =  ACOgIP) 

CALL  ELOWMI  \JI  fits.  Awes,  SRC  ,  SNK,  I  ,  J,HI  ,r  l  ffc  ,NC,NF,NA,NH,CHJ  ) 

DO  3  Jl=  I,  ARCS 
P  =  ICHGIMU  J  I  )  I 
HI!  Jl  )=  L 
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P  =  BCHGt FLOW! Jill 
)  FLOW  I J I  )=  L 

CO  4  11=  1, NODES 
P  =  PCHGINFI  I  I  I  I 
MFC  I  I  I  =  L 
P*  BCHGINCI  HI) 

<.  NCI  ID*  L 
P  =  HCHGIORJ) 

oej  =  l 

RETURN 

END 

HHF1C  FLOWM 

SUP ROUT  I  NT  FLOWN  < NOUE $ , ARC S  ,  SRC  ,  SNK , I , J , H I , F LOW, NC , NF , N A , NH , OH J ) 
INTI Gl R  NOONS, ARCS, SNK, SRC, I  I  2  COO ) , J I  2000  I ,H I ( 2000  I , F LOW  I 2G00 l 
INTEGER  NCI  1  COO  I , NF ( ICOC) , NA ( 1 000  I . NBI 1000 ) , OB J 
C 
C 

C  DEFINITION  OF  CALLING  SFCUENCE 

C 

C 

C  NAME  USF 

C 

C  NODES  NUMBER  OF  NODES 

C  ARCS  NUMBER  OF  A*CS 

C  SRC  SOURCE  NODE. 

C  SNK  SINK  N(jn 

C  I  LIST  (IT  FIRS!  NOOLS 

C  J  LIST  OF  StCONO  NODES 

C  HI  UPPER  HOUND  S  FDR  ARCS 

C  FLOW  AMOUNT  OF  FLOW  IN  ARCS 
C  NC  NODE  CAPACITY  I  INPUT  I 

C  NF  NODE  FLOW 

C  NA  SOURCF  NODE  LABELS 

C  NH  SINK  NODC  LABELS 

C  OH  J  OBJECTIVE  VALUE 

C 

L  BEGIN 

INTEGER  A »  A  A  » N ,  N 1 » N2 » G  T  ARC S, INC,  LABEL 
LOGICAL  TYPF 
GTARCS  =  ARCS  »l 

no  m  a  *  l , arc  s 

FIOWIAI  =  0 

ir  I  1(A). LF.O. OR. 1(A). GT. NODES)  GO  TO  999 
IF  I JIA). Lfc.C. OR. JIAI.GT. NODES)  GO  TO  999 
IF  I IIAl.Fg.JlAI)  GO  TO  999 
HI  I  A)  =  I  A  0  S I H I  (  A  I  ) 

10  CONTINUE 

DO  ?C  N  =  1,  NODES 
NFINI  =  0 

NCINI  =  lABSINCINl) 

2C  CONT INUt 
OBJ  =  G 

C  /< RO  NODE  L  ABELS 

100  LP  200  N  =1, NODES 
NA(N)  =0 
NHINI  =0 
2 OP  CONTINUE 
C  I  API  l  SOURCE  NODI 
KOI  SRC)  =  1 

IE  (NCI  SRC  I.FC.NF  I  SRC)  )  RETURN 
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NAISRCI  =  GT  ARCS 
C  LABEL 

210  LABEL  =  0 

00  250  A=l,ARCS 

N 1  =  1(A) 

IE  (  N 1  •  L  T  .  0 )  (.<:  Tf)  250 
N2  =  J (  A  ) 

IF  (FLOW! A ).L I .0)  GO  TO  225 
IF  (  N  A  (  N  l  )  .NF.OI  GO  TCJ  220 
IF  (NBIN2) .LG. 01  GO  TO  250 
IF  (FLOMA  I.EU.O)  GO  TO  224 
N  A  I  N 1  I  =  -A 

IF  (NR(Nl).NE.'J)  GO  TO  240 
NR(Nl)  =  -GTAKCS 
GO  10  24C 

220  IT  CNRIN2) .Nf .O.OR.FLCWl A) .FO.HI (A) )  GO  IC  245 
NR ( N  2 )  =  A 

IF  (NF(N2).E0.NC(N?) )  GO  TO  240 
NAIN2)  =  GIAKLS 

GO  TO  24U 

224  IF  (NAIN2) .r^.O)  GO  TO  250 
GO  TO  226 

225  IF  (NAIN2l.f-Q.OI  GO  TO  230 

226  IE  (NR(Nl).NE.O.OR.FLOW(AI.EO.(-Hl (A)))  GC  10  245 
NO  (  N  I  )  =  -A 

IF  (NF(Nl).FC.NC(Nt)  )  GO  TO  240 
NAINl)  =  ♦GTAHCS 

GO  TO  240 

230  IF  (NRINl ) .Ft. O)  GO  TO  250 
NA(N2)  =  ♦ A 

IF  INDIN?) .NE • 0 1  GO  TO  240 
NP ( N? )  =  -GTAKCS 
240  L ARE  L  =  l 

IF  ( NA I  SNK  I.NF.  .0)  GO  TO  260 
245  1(A)  =  -N l 

250  CONTINUE 

C  GO  BACK  AND  LABEL  NORF  IF  SOHE  NODE  WAS  LARELEC  ON  L  A  s  T  LOOP 
IF  (LABEL. NF.C)  GO  TO  210 
C  RESTORE  POSITIVE  SIGNS  TO  FIRST  NODF  LIST 
260  00  270  A  =  l, ARC  S 

1(A)  =  I  AH  S (  1(A)) 

270  CONTINUE 

C  IE  NOTHING  LABELED  ON  LAST  LOOP,  DONE 
IF  (IARFL.EU.CJ  RETURN 
C  BREAKTHRU.  FIND  THE  INCREMENT 
300  INC  «  NCI  SMC I-NF ( SRC  ) 

C  FOLLOW  PATH  RACK  FROM  SINK 
N  *  SNK 
TVPF  =  .  TRUE  • 

310  A A*  NA ( N ) 

IF  (  .NO  I . TYPE )  A  A  =  NR  (  N  ) 

A  r  IARS(AA) 

IF  (A.GT.AKCSI  GO  TO  320 
IF  (AA.LT.O)  GO  TO  315 
N2  =  1(A) 

IF  (FLOW! AI.LT.O)  GO  TO  316 
INC  =  MINOl  INC, HI  (AI-FLOW(A)  ) 

GO  I 0  318 
315  N2  =  J ( A ) 

IF  (FLOW! AI.Lf .0)  GO  TO  317 


t 
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Me 

INC 

- 

HMtM 

1  NO , 

1 AHSI TLOWI A  >  1  ) 

no 

Til 

11 '1 

■»!  / 

INF. 

r 

M  I  ",  0  ( 

IM  , 

II  Owl  A  )  ♦  )<  I  (All 

J  IB 

N  - 

N  2 

GC 

rr 

3  50 

II 

(TYPE  i  <; 

i:  T i 

1^5 

INC 

MIN01 INC, 

•  1  ll,M 

GO 

IN 

350 

3?5 

INC 

- 

M  INOI 

1  NC  , 

NC IN) -NF (Nil 

350 

TYPE  - 

.NOT 

.  I  YPI 

1  F 

<r. . 

M  .  SR 

0  )  0 

1)  1C  310 

f.  lUCKfcMIN!  ARCS 


V  =  SNK 

I  YPl  =  .  TKUI  . 

UMJ  =  UliJ  ♦  UX 
ISC  A  A  =  NO  IN) 

IF  (  TYPE  )  AA  --  *  A  I  N  ) 

A  =  I  A  H  S  (  AA  ) 

If  (A. NT. ARCS)  (.»)  TO  36S 
If  (AA.LI.U)  Gli  IN  355 
FlllWlAJ  =  fli'MAI  ♦  INC 
■Y  =  1  (  A  I 

(  c  in  3 7u 

365  nilUAI  -  Ki’klAI  -  INC 
W  =  J  (  A  ) 
r,c  in  3  7o 

365  IF  (  TYPF  )  (,il  IH  36* 

NFINl  =  NFI  J»  -  INC 

GG  Id  3  7c 

16  H  NF(N)  =  M  F  C  N I  ♦  INC 
3/0  TYPI  =  .  Nil  I  .  T  YPfc 

If  I  \.\IF.SRC.nR.  IYPE  )  GO  1G  350 
C  FUJI*  iNCiUMtNm),  REllKN  TO  LABELING 
(.0  in  100 

05^  PR  1  .1  NNfi,  A,  I  (  A)  ,  J(  A  I 

NOM  (  OKI'  A  I  (  5dH  ARCS  INCOHRCCILY  SETUP  IN  NCTwCRK  ROUTINE,  ARC  |6,3H 
> 16, 1H, 16, I H ) 1 
KFTIHN 
ENG 
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Appendix  B 
EXAMPLE 

For  illustrative  purposes,  the  algorithm  is  used  to  find  the 
optimal  placement  of  forces  for  the  network  of  Fig.  1,  with  output 
printed  after  each  iteration.  Ten  forces  are  available,  and  all  arcs 
and  nodes  have  identical  probabilities  associated  with  them.  Thus, 
for  each  node,  i,  and  each  arc  (i,j),  we  have: 


p(i)L 

S 

.  666b 

PU»  J)  j 

= 

.6666 

p(i)2 

= 

.8332 

P(i.  j>2 

= 

.8332 

P(i)  3 

= 

.8748 

P(  i *  j ) 3 

= 

oc 

'■'■J 

oc 

P<1>4 

= 

.8852 

S 

.8852 

P(i>5 

= 

.8878 

pU.j)5 

- 

.8878 

P(i)6 

= 

.8884 

p(i .  j)6 

= 

.8884 

P(i)? 

= 

.8886 

PU,  j>7 

= 

.8886 

T3 

'w' 

00 

B 

.8887 

pd,j)8 

= 

.8887 

P(i)9 

= 

.8888 

pU,  j)9 

E 

.  8b88 

p<1,10 

= 

.8889 

= 

.8889 

The  input  deck  is  given  in  Fig.  2  and  the  output  deck  or  results 
follow . 


-28- 


0090150010090101 

00 l 00 l UO l 002002 00 3 00 3004 005005 

00600o00700700f 

00200300400300500400500 i006008 
00 7009008009009 

0 . 66660 .83320. 874H0. 80520.88780. 88840. HRB60.B8870. 88800. 8889 
0.66660. 8  3  320. 8  7480. 88  *320.  88  780. 8884  0. 88  860. 808  70. 88880. 888  9 
0.66660.83320.87480.88520.88780.88840.88860.88870.88880.0889 
l -.66660. 8 3 320. 8  74  80. 88  520. 88  780.88840.88860.8  88  70.88880.8889 
O. 66660. 8 3320. 8 7480. 88 520. 88 780. 88840. M806O. 888 70. 88800. 888 9 
O.6666O.83320.97480.88520. 88780 . RR840 . 88 860 . 888 70 .88880 . 8889 
0. 66660. 8 3 320. 87480. 88 520. 88780. 88840. H8 860. 888 70.88880. 8889 
O.6666O.83320.8748O.88520. 88 78C . 88840 . 88 860 . 88 8 70 .88880 . 8889 
0.66660.8  3  320.8  74 00. 88520. 08 780. 88840.80 860. 888 70. 88 880. 8889 
0.66660. 8 3 320. 8 7480. 88520. 08 780. 88840. 88 860. 88 8 70. 88880.8889 
0.66660.8  3  320. 0  74 80. bB 5 20. 88  780. 88840.88 86 0.888 70. 88880. 8889 
0.66660. 8 3 320. 8 7480. 88 520. 88 780. 88840. 88 860. 888 70. 88880. 8889 
0.66660.83320.87480.88520. 88780.88840.88860.88870.88880. 8889 
o. 66660. 8 3 320. 8 74 00. 88 520. 88 7 80. 88 84 0.80 8 60. 888 70. 88 880. 8889 
0.66660.8  3  320.8  748  0. 80  520. 88  780 . B8840 . HH860 . 888 70 . R8B80 . 88«o 
0. 66660. 8 3 320. 8 74  ho. 8ft 520. 88 780. 88840. P0 860. 8 88 70. 88 880. 8889 
<  . 66660. 8 3 320. H 74 8 0.88 5 20. 88 7 80. 88 840. 88 860.888 70. 88 880. 8809 
0.66660.8  3 320. H 74 8 O. 88 5 20. 88 780. 88 840. H8 8 60. 8RP7O.R0 880. 8809 
0.66660.8  3  320.  H  748  0. 8H  5  20.  08  780. 88  84  0  . »'  8  86  0 . 088  70 . 888H0 . 8889 
(<.  66660. 8  3  320.  H  748  O.HH  5  20.  88  700. 88  040 . 88  860  .  R88  70  .  RBBRO .  8889 
0.66660.8  3  320.8  7480.  ft H 520 . 88 780 . 88840 . 8H 860 . BOH  70 . 8B8H0 . 8889 
0.66660.8  3  320. H  7480.88  520. 8R 780. 0R 340. 8  8 860. 8 88 70.88880. R8R9 
0 .66660.83320.87480. 8H 520. 88 780. Bft840.HH 860. 88P70. 88880. 8889 
6.<  6660. 8  3  320. 8  74 HO.  88  520.  8  8  780 .  HR 84  0  .  l'H 860 . 88  8 70 . 8ft 880 . 880  9 
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Fig.  2--Input  deck 


OUTPUT  DECK 


SOURCE  NTT  FLOW  SINK 


1  .*>0000 


MAX.  PATH  VALUE 


0.  3333 


FLOW 


CUT  SFT 


FORCES 


1  .*>0000 
1.50000 
o.r-oooo 
0.00000 
1  .*>('000 
1.50000 
0.00000 
O.OuUOO 
1  .>0000 


l.OOCOO 

0.00000 

O.OOGCO 

o.oooco 

o.oooco 

o.oooco 

0.00000 

0.00000 

0.00000 


Fl.uW 


CUT  SET 


FORCES 


3 

5 


5,  H 

6,  1 
h.  9 
7,  H 


1 .50000 
0.00000 
0.00000 
0.00000 
1.50000 
0.00000 
0.00000 
0.00000 
l  .50000 
0.00000 
0.00000 
1.50000 
0.00000 
0.0(000 
0 .00000 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.oooco 

0.00000 

0.00000 

0.00000 

o.oocoo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 
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S'UKCF  r. f  T  FLIJW  SINK 

l  I  . snooo  q 


in 

FI  Ow 

f.u  1  S  t  T 

i 

1  .'>0000 

0 

s 

l  .‘>0000 

u 

< 

0.00000 

0 

4 

0.00000 

0 

s 

1  .SOOOti 

l 

i 

l  .SOOOG 

0 

7 

0.00000 

0 

i 

o.orooo 

0 

') 

1  .‘>0000 

0 

4 

f 

FLOW 

CUT  SLT 

1 . 

1  ,'OlUOv 

0 

1 . 

< 

0.00000 

0 

1 . 

4 

0.00000 

0 

•  '  » 

) 

o.oooou 

u 

> 

'  • 

s 

1  .  >0000 

0 

5  t 

4 

o.oooou 

0 

0.00000 

0 

'< . 

> 

0.00000 

0 

'*  * 

.*> 

1  .50000 

0 

t 

It 

C.UUOOO 

0 

o  » 

7 

O.UuOOO 

0 

'  • 

) 

1  .SPOOO 

0 

< , 

vt 

0.00000 

0 

/ , 

') 

0.00000 

0 

% 

0.00000 

0 

MAX.  PATH  VALUF 
O.lllll 


FUHLfS 


i .cncoo 
o.occcu 
0.00000 
O.OCCCU 

l .occoo 

O.COCCO 

o.occco 

O.OPCOO 

o.occco 


FOKCI  S 


O.OCOOO  I 

o.occco 

o.uocoo 

o.occoo 

o.cccoo  i 

o.occco 

o.cccoo 

o.occco 

o.orooo  \ 

o.occco 

o.coooo  ! 

o.occco 

0.00000 

o.ococo 

0.00000 
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SmiKC.f  Nf  T  FLOW  SINK 


MAX. 


l  i.sonon  y 


NODE  FLOW 


cur  SET  FORCES 


1 

? 

3 

A 

‘j 

0 

7 

h 

■V 


1  .-jOOOO 
l . >0000 
0.00000 
0.00000 
1 • SOOOO 
l .sooon 
0.0000( 
0.00000 
1  .  >  U  0  (j  u 


0 

0 

0 

0 

0 

0 

0 

0 

L 


1.00000 
o.oocoo 
o.occoo 
o.ccooo 
1 .oocco 
o.occcr 
0  .  UCCOC 

0 . orcoo 
l .( ooco 


ARC 


FLOW 


(UT  SET  FORCES 


1 . 

/ 

1 . >0000 

1 . 

5 

o .uuuoo 

1 . 

A 

0.00000 

f . 

3 

U  .('0000 

/ . 

T 

1  .‘>0000 

3. 

‘i 

0.00000 

3, 

s 

0.00000 

A  t 

> 

0.00000 

’*  * 

<s 

1  .  >0000 

> « 

8 

0.00000 

7 

o.uuooo 

6, 

7 

1 .SOOOO 

7, 

8 

0.00000 

7, 

4 

0.00000 

) 

0.00000 

0 

0 

o 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0. oocco 
0. oocco 
0.00000 
0.0001  c 
0.00000 
o.oocoo 

O.OOCOO 
0. OCOCt 
0.00000 
0.00000 
O.OCOCO 
0.00000 
0.00000 
0.00000 
0.00000 


PATH  VALUE 
0.03704 
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SOURCF  NET  FLOW  SINK 
1  2.00098  9 


J'MJF 

Flow 

C'JT  SET 

1 

2.00098 

1 

2 

1  .50000 

0 

3 

0. 5009ft 

0 

A 

0.00000 

0 

5 

2.0009H 

0 

6 

1.50000 

0 

7 

0.5G09H 

0 

q 

0.50098 

0 

9 

2.00098 

0 

A9C 

FKJW 

CUT  SET 

1  . 

2 

1 .50000 

0 

1  . 

3 

0.50098 

0 

1  , 

A 

0.00000 

0 

2  , 

0.00000 

0 

2. 

5 

1.50000 

0 

3. 

A 

0.00000 

0 

3. 

5 

0.50090 

0 

A  t 

) 

0.00000 

0 

5  * 

G 

1 .50000 

0 

>  * 

8 

0.50098 

0 

t>  • 

7 

0.00000 

0 

9  » 

9 

1 .50000 

0 

q , 

7 

0.50098 

0 

r . 

lV 

0.50098 

0 

H  , 

V 

0.00000 

0 

MAX.  PATH  VALUE 
0.01653 


FOWL  FS 


2.00000 
0.00000 
0.00000 
0.00000 
1 .00000 
0.00000 
0.00000 
0.00000 
1 .00000 


FORCES 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 
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SOURCE 

NFT  FLOW 

SINK 

l 

2.00090 

9 

NC10E 

FLOW 

CUT  SfcT 

1 

2.00098 

0 

2 

l  .80000 

0 

3 

0.80098 

0 

4 

0.00000 

0 

8 

2.00098 

1 

6 

1  .8  0000 

0 

7 

0.80098 

0 

8 

0.80098 

0 

9 

2.00098 

0 

ARC 

F  LOW 

CUT  SET 

1  . 

2 

l .80000 

0 

1. 

3 

0.80098 

0 

1  . 

4 

0.00000 

0 

2  t 

3 

0.00000 

0 

?. 

8 

I .80000 

0 

1, 

4 

0.00000 

0 

3, 

8 

0.80098 

0 

9  t 

8 

0.00000 

0 

8, 

6 

l .80000 

0 

8, 

8 

0.8009R 

0 

t 

7 

0.00000 

0 

6, 

9 

1.80000 

0 

H, 

7 

0.80098 

0 

7, 

9 

0.80098 

0 

«» 

9 

0.00000 

0 

MAX.  PATH  VALUE 
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Appendlx  C 

CONVEXITY  ASSUMPTION  ON  MULTIPLE  FORCES 


It  was  mentioned  earlier  that  the  assumption  pC i-) ic-t-I  "  P^k  S 
p(i)^  -  p( i ) ^  ^  assures  that  allocating  [rr(i)]  +  1  forces  with  probability 
ff(i)  and  [n(i)]  forces  with  probability  1  -  tt ( i )  minimizes  the  proba¬ 
bility  that  the  infiltrator  can  successfully  cross  node  i  for  given 


n(i) .  This  appendix  will  verify  that  statement. 

Let  P[k,l]  be  the  probability  that  k  forces  are  placed  at  node  i. 
Then,  EkP{k,i}  =  rr(i)  and  the  probability  associated  with  node  i  is 
1  -  ZP{k,i}p(i)^.  Minimizing  this  probability  is  equivalent  to 
maximizing  EP[k,i}p(i)  .  The  following  theorem  shows  that  for  given 
n(i),  this  occurs  when  the  conditions  above  are  satisfied. 

THEOREM.  The  quantity  (i)EP{k, i] p(i)^  is  maximized,  subject  to 
the  constraint  (ii)£kP{k, i}  =  rr(i)  when  P{[n(i)]  +  l,i|  =  ft(i), 
P([n(i)],i}  =  1  -  fr(i),  and  all  other  P{k,i}  =  0. 

PROOF .  Let  h  =  max(k/P{k,i}  >  0} ,  i  =  min{k/P{k,i}  >  Oj,  and 
d  *  h  -  l.  Let  P[k,i)  =  P{k,i}  be  a  set  of  values  of  smallest  d  which 


maximizes  (i)  subject  to  (ii)  .  Suppose  the  theorem  is  false.  Then 
h  -  £  =  d  £  2.  Let  m  =  min{P[h,i},  P( £ , i} J  .  Set  P[h,i}  =  P{ h, i}  -  m 
and  P[h  -  1 , i }  =  P{h  -  1 , i }  +  m.  Then  set  P[ji,i}  =  P{£,ij  -  m  and 
increase  P [i  +  l.i}  by  m.  (If  h  -  l  =  2,  then  PU  +  1  =  ?{t  +  1}  +  2m). 
Either  p£ h , i }  or  P{£,i}  is  now  zero,  decreasing  d  by  at  least  one  unit. 
The  increase  in  (i)  is  equal  to  m(-p(i)2  +  p(Oj1+1  -  p(i>h  +  p(i)h_  ). 
Since  (ii)  is  still  satisfied,  this  is  strictly  negative  from  our  choice 
of  P[k,t}.  However,  since  h  >  l  +  1, 
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P(i)h  -  P<i>h.i  *  P<1>£+1  ■  P<i)£ 

-  P(i)£  +  p(i)m  -  P(i)h  +  PU),,.!  -*  0 
for  a  contradiction.  QED. 
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Appendix  D 

FINDING  MAXIMUM  PATH  VALUE 


Thi3  appendix  briefly  explains  the  procedure  for  keeping  track 
of  the  maximum  path  value  during  the  algorithm. 

Let  K  be  the  maximum  path  value  at  the  beginning  of  the  algorithm 
and  let  V  be  the  value  of  C,  the  minimum  cut  set  .  Substituting  the 
generalized  capacities  (Expression  (4a-b))  for  the  ones  (Expression 
(2a-b))  that  appear  in  Expression  (3),  one  obtains: 

E  An(_l)  ±  E  An(i  ._i)  _  i  v 

AK  ~  K 


But  £  An(i)  +  £  Arr(i,j)  *  MV,  where  M  is  as  defined  in  step  5  of  the 
multiple  intercepting  force  algorithm.  Thus 


MV 

AK 


or  K  -  AK  -  K( 1  -  M) . 
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